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Abstract 

We study approximation of embeddings between finite dimensional 
L p spaces in the quantum model of computation. For the quantum 
query complexity of this problem matching (up to logarithmic fac- 
tors) upper and lower bounds are obtained. The results show that for 
certain regions of the parameter domain quantum computation can 
essentially improve the rate of convergence of classical deterministic 
or randomized approximation, while there are other regions where the 
best possible rates coincide for all three settings. These results serve as 
a crucial building block for analyzing approximation in function spaces 
in a subsequent paper 



1 Introduction 

In this paper we continue the investigation of numerical problems of analysis 
in the quantum model of computation. In a number of papers the integration 
problem and its discretized version, the mean computation, were studied 
and matching upper and lower bounds (often up to logarithmic factors) 
were established for various function classes. See the references [21 EJ EE 

imiElimiiniESIISIE]. lt turned out that for these tyP es of problems 
quantum computing can reach an exponential speedup over deterministic 
classical computation and a quadratic speedup over randomized classical 
computation. 

All these problems are such that the solution is a single number. There- 
fore the question arises what happens if we consider problems whose solu- 



1 



tion is a family of numbers or, in other words, a function. A particularly 
typical situation is function approximation, where the solution is just the 
input function itself (and we are asked to compute an approximation to it 
in a given norm). A first consideration of an approximation problem in the 
quantum setting appears in [20], but no matching upper and lower bounds 
were obtained. 

In the present paper we provide the first results for approximation in the 
quantum model of computation, with matching upper and lower bounds. We 
start with the very basic situation: the approximation of the embedding J^ q 
of Lp into Lq, or, in other words, the approximation of A-sequences with 
bounded Lp norm in the norm of L q . These embeddings are the elementary 
building blocks of embeddings of function spaces - in the same way as mean 
computation is the elementary building block of integration, see |1(J| and 
for more on this principle. Our results show that for p < q, the quantum 
model of computation can bring an acceleration up to a factor A -1 of the 
rate of the classical (deterministic or randomized) setting. On the other 
hand, for p > q, the optimal rate is the same for all three settings, so in 
these cases there is no speedup of the rate by quantum computation. 

We prove that the following version of Grover's quantum search algo- 
rithm is optimal: we find all coordinates of / € Lp with absolute value not 
smaller than a suitably chosen threshold and set the other coordinates to 
zero. The crucial new element in proving lower bounds is a multiplicative 
inequality for the n-th minimal query error, which is - in a wide sense - 
analogous to multiplicativity properties of s-numbers, see [2T] . 

In a subsequent paper we show that, similarly to the analysis in 
[TU] . sufficiently precise knowledge about the embeddings Jp q leads to a 
full understanding of the infinite dimensional problem of approximation of 
functions from Sobolev spaces. 

The paper is organized as follows. In Section 2 we recall notation from 
the quantum setting of information-based complexity theory as developed 
in [Hj. In Section 3 we derive some new general results which will be needed 
later on. Section 4 contains the main results on approximation of embed- 
dings of Lp into Lq spaces. Finally, in Section 5 we give some comments 
on the quantum bit model and a summary including comparisons to the 
respective results in the classical deterministic and randomized setting. 

For more details on the quantum setting of information-based complexity 
we refer to [HJ, to the survey ^H] , and to an introduction . For the classical 
settings of information-based complexity theory see |181 1261 |7]. General 
background on quantum computing can be found in the surveys HJ [2l] 
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and in the monographs [221 E] and . 



2 Notation 

For nonempty sets D and K, we denote by J-{D, K) the set of all functions 
from D to K. For a normed space G we let 0(G) = {g G G| ||g||c7 < 1} 
denote the unit ball of G. Let .F C J-(D, K) be a nonempty subset. Let 
K stand for either R or C, the field of real or complex numbers, let G be 
a normed space over K, and let S : F — > G be a mapping. We seek to 
approximate <S(/) for / G F by means of quantum computations. Let -Hi 
be the two-dimensional complex Hilbert space C 2 , with its unit vector basis 
{e ,ei}, let 

H m = Hi® •••<g> ill, 
s » ' 

m 

equipped with the tensor Hilbert space structure. Denote 

Z[0,iV) := {0,...,iV-l} 

for G N, where we agree to write, as usual, N = {1,2, ... } and No = 
N U {0}. Let C m = {\i) : i G Z[0, 2" 1 )} be the canonical basis of H m , where 
\i) stands for e JO <S> ■ ■ ■ ® ej m _ 1 , z = YlT^o 3k2 m ~ 1 ~ k is the binary expansion 
of i. Let U{H m ) denote the set of unitary operators on H m . 
A quantum query on F is given by a tuple 

Q = (m, m! , m", Z, r, /?), 

where m,m' ,m" G N,m' + m" < m, Z C Z[0, 2 m ) is a nonempty subset, 
and 

r : Z -> D 
/3 : A'^ Z[0,2 m ") 

are arbitrary mappings. We let m(Q) := m be the number of qubits of Q. 

Given a query Q, we define for each / G F the unitary operator Qf by 
setting for |i) \x) \y) G C m = C m i 8) C m // ® C„ 



"m—m —m" 



Qf 



l) IX) 



\i)\x®(3(f(r(i))))\y) if < G Z 
|i) |x) |y) otherwise, 



where © means addition modulo 2 m . 

A quantum algorithm on F with no measurement is a tuple 

A=(Q,(Uj)] =0 ). 
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Here Q is a quantum query on F, n G No and C/j G U{H m ) (J = 0, . . . , n), 
with m = m(Q). Given f £ F, we define ^4/ G lA{H m ) as 

A f = U n Q f U n - 1 ...U 1 QfU . 

We denote by n g (j4) := n the number of queries and by m{A) = m = 
m(Q) the number of qubits of A. Let (Af(x, y)) x ,y ez\o,2 m ) be the matrix 
of the transformation Af in the canonical basis C m , that is, Af(x,y) = 
(A f \y),\x)). 

A quantum algorithm from F to G with k measurements is a tuple 

A=((A l )^,(b l )^,<p), 

where k G N, Ai (/ = 0, . . . , k — 1) are quantum algorithms on i 7 with no 
measurement, 

6 GZ[0,2 m »), 

i-i 

6, : JJZ[0,2 m< ) -> Z[0,2 m <) (1 < Z < fc - 1), 

i=0 

where m; := m(j4j), and 

fc-i 
z=o 

The output of ^4 at input / G -F is a probability measure -A(/) on G, defined 
as follows. First put 

PAj(x ,... = |A 0)/ (a;o,6o)| 2 |-4i,/(^i,fei(a;o))| 2 ■ • • 

. . . |vl fe _ij(x fc _i A_i(x , • • • ,x k _ 2 ))\ 2 - 

Then define A(f) by setting for any subset CCG 

A{f){C) = ^2 p A ,f(x ,...,x k -i). 

ip(x ,-,x k -i)ec 

Let n q {A) := Yl*i=o n q(Ai) denote the number of queries used by A. For 
more details and background see jSj. Note that we often use the term 
'quantum algorithm' (or just 'algorithm'), meaning a quantum algorithm 
with measurement (s). 

If A is an algorithm with one measurement, the above definition simpli- 
fies essentially. Such an algorithm is given by 

A = (A , bo, <p), A = (Q,(Uj)] =0 ). (1) 
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The quantum computation is carried out on m := m(Q) qubits. For f £ F 
the algorithm starts in the state |&o) an d produces 

\il>f) = U n Q f U n - 1 ...U 1 Q f U \b ). (2) 

Let 

2 m -l 
i=0 

(referring to the notation above, we have cii / = Ao /(i, bo))- Then A outputs 
the element ip(i) G G with probability |a,j| 2 . It is shown in [Sj, Lemma 1, 
that for each algorithm A with k measurements there is an algorithm A with 
one measurement such that A(/) = A(f) for all / G F and A uses just twice 
the number of queries of A, that is, n q (A) = 2n q (A). Hence, as long as we 
are concerned with studying the minimal query error (see below) up to the 
order, that is, up to constant factors, we can restrict ourselves to algorithms 
with one measurement. 

Let 6 > 0. For a quantum algorithm A we define the (probabilistic) error 
at / G F as follows. Let ( be a random variable with distribution A(f). 
Then 

e(S,AJ,9) = inf {e > | P{||5(/) - CII > e} < 9} 

(note that this infimum is always attained). Hence e(S,A,f,9) < e iff the 
probability that the algorithm A computes S(f) with error at most e is at 
least 1 — 9. Observe that for algorithms with one measurement (©-©), 

p{\\s(f) - C|| > e} = y, K/i 2 - 

i: \\S(f)-<p(i)\\>e 

Trivially, e(5, A, f,9)=0 for 6 > 1. We put 

e(S,A,F,^) = supe(S,A,/,0) 

(we allow the value +co for this quantity). Furthermore, we set 

= inf{e(S, A, F,9) \ A is any quantum algorithm with n q {A) < n}. 

It is customary to consider these quantities at a fixed error probability level. 
We denote 

e(5,A,/)=e(5,A/,l/4) 
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and similarly, 

e(S,A,F) = e(S,A,F,l/4), e^(S, F) = e^(S, F, 1/4). 

The choice 9 = 1/4 is arbitrary - any fixed 9 < 1/2 would do. The quantity 
en(S, F) is central for our study - it is the n-th minimal query error, that 
is, the smallest error which can be reached using at most n queries. Note 
that it essentially suffices to study en(S, F) instead of en(S, F, 9), since with 
0{u) repetitions, the error probability can be reduced to 2~ v (see Lemmas 
13 H and Corollary □below). 

3 Some General Results 

Let G and G be normed spaces. Recall that a mapping <I> : G — > G is said 
to be Lipschitz, if there is a constant c > such that 

||3>(a;) — <&(y)||g < c || x — y\\a for all x,y G G. 

The Lipschitz constant H^HLip is the smallest constant c such that the rela- 
tion above holds. 

Given a quantum algorithm A from F to G and a mapping <I> from G to 
G, the algorithm <E> o A is defined as the composition, meaning that ip in the 
definition of A is replaced by $ o ip (this is a special case of the definition 
of the composition given in jSj, p. 13). The following direct consequence of 
the definitions will be needed later. 

Lemma 1. Let S be a mapping and A a quantum algorithm, both from F 
to G. Let & be a Lipschitz mapping from G to G. Then for each f G F and 
9>0, 

e($oS,*oA,/,0) < \\nu P e(S,A,f,9). 

Consequently, for each n € N, 

eU^oS,F,9)<\M Up e^S,F,9). 

The next result was shown in jlUj . Lemma 2, for G = R, but the proof 
of the general case is identical to that one. 

Lemma 2. Let D,K and F C F(D,K) be nonempty sets, G a normed 
space, let k G No and let Si : F — > G (I = 0, . . . ,k) be mappings. Define 
S:F^Gby S(f) = Ef =0 Si(f) (f€F). Let , . . . , 9 k > 0, n , . . . , n k € 
No and put n = Ya=o n i- Then 

k k 

el{S,F^i) <Y.<^ F ^)- 

1=0 1=0 
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The following results are generalizations of the usual procedure of "boost- 
ing the success probability", which decreases the failure probability by re- 
peating the algorithm a number of times and computing the median of the 
outputs (see, e.g., jS], Lemma 3). This works for algorithms whose outputs 
are real numbers. Since there is no natural linear order on a normed space, 
in general, the latter step has to be changed suitably when dealing with 
outputs in a normed space G. 

For this purpose, let v £ N. Let /j, : R 1 ' — > R be the mapping given by 
the median, that is, /i(ao, . . . ,a u -\) is the value of the [~(^ + l)/2]-th element 
of the non-decreasing rearrangement of (aj). First we deal with the case that 
G is a space of the form G = Zoo CO > where T is a nonempty set and Zoo CO 
denotes the space of all bounded real-valued functions on T, equipped with 
the supremum norm H^H^m = su Ptgr Define p, : l 00 (T) v — > Zoo CO 

as follows: 



that is, we apply the median componentwise. For any algorithm A from 
F to Zoo(T) denote by Jjl{A v ) := p(A, ...,A) the composed algorithm (see 
again p. 13 of of repeating v times the algorithm A and applying p, to 
the outputs. 

Lemma 3. Let A be any quantum algorithm and S be any mapping, both 
from F to l^T), and let u EN, Then for each f £ F, 



Proof. Fix / £ F. Let £q, . . . , C^-i be independent random variables with 
distribution A(f). Let \i be the indicator function of the set 



Then P{xi = 1} < 1/4. Hoeffding's inequality, see, e.g., p. 191, yields 



n{gQ,...,g v -i) = (fi(g (t),...,g u -i(t)))teT 



e(S,p,(A»)J,e 



'/ s )<e(S,A,f). 



{\\S(f)-Ci\\u(T)>e(S,A,f)}. 



K i=0 ) \ i=0 ) 

Hence, with probability at least 1 — e~ u ^ 8 , 

\\ s (f) ~ Cilhoo(T) < e(S,A,f)}\ > v/2. 
It follows from © that for alH G T 




(4) 



\{i | \S(f)(t)-Ci(t)\<e(S,AJ)}\>u/2. 
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Consequently, 

W)(t) - KCo(t), (u-i(t))\ < e(S, A, /), 
which means that 

||5(/)-A(Co,... ! C,-i)lb 00 (r ) <e(5,A/)- (5) 

Since ® holds with probability at least 1 — e~ u ^ 8 , so does ©. □ 

Now let G be a general normed space. For the following construction 
we consider G as a space over R (each normed space over C can also be 
considered as a normed space over R). We define for each 5 > a suitable 
mapping tpg : G v — ► G as follows. Let G* denote the dual of G, that is, the 
space of all bounded linear functionals on G. Let T C B{G*) be a norming 
set, i.e., for all g £ G 

\\g\\ = sup (6) 

(such a T always exists, one can take, for example, T = B(G*) itself). 
Then G can be identified with a subspace of loo{T) via the embedding map 
J : G — > loo{T) defined by 

Jg = (t{g))teT 

in such a way that the norm is preserved: || JgWi^^) = IIs'IIg- For 5 > 0, let 
finally its : loo{T) — > G be any (^-approximate metric projection, by which 
we mean a mapping satisfying 

\\ x - ^OOlLcr) < (i + ^mf lk-5'lb 00 (r) 

for all x G loo{T). We define ^ : G u G by setting 

1p 5 = TT S flo J U . 

Lemma 4. Lei A be any quantum algorithm and S be any mapping, both 
from F to a normed space G, let v £ N and 5 > 0. Then for each f G F, 

e(S, MA U )J, e~ u/8 ) < (2 + S)e(S, A, /). 

Proof. It follows from © that 

e(J5, JA/) = e(5,A,/)- 



S 



By Lemma 01 

e(JS, fi((JA) v ), f, e~ v ' 8 ) < e(JS, J A, f) = e(S, A, /). 

Let C be a random variable with values in Z 00 (7") with distribution /x(( J A) u )(f). 
Then with probability at least 1 — e~ u ^ s , 

\\JS(f)-(\\<e(S,A,f). 

Hence 

\\s(f) -MOW < IW)-CII + IIC-^(C)[| 

< e(S,A,f) + (l + 5)\\C-S(f)\\ 

< (2 + 5)e(S,AJ). 

But tts(C) is a random variable with distribution 

7T S of,((JAY)(f)=MA u )(f), 

and it follows that 

e(S,MA u )J,e~ u/8 ) < (2 + 5)e(S,A,f). 

□ 

Corollary 1. Let S be any mapping from F C T(D,K) to a normed space 
G. Then for each n, v G N, 

eUS,F,e-»/ 8 )<2eUS,F). 

If G is a space of the form Zoo CO f or some setT, then the constant 2 above 
can be replaced by 1. 

The definition of ip$ is not constructive (and neither is that of p,, if T 
is infinite). Since we are dealing with the quantum query complexity, the 
cost of (classically) computing the (p part of a quantum algorithm (this is 
the place where jl and ip$ enter) are generally neglected. However, if one 
looks for a more efficient procedure, here is one which leads to the constant 
3 instead of 2 + 5 of LemmaHJ Define g : G v — > G as follows: g(go, ■ ■ ■ , 9v~\) 
is the element g^, where 

i = argmin/i([|3o-0i[|)---,||&/-i ~ 
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(if there is more than one index % at which the minimum is attained, we 
choose the smallest index, just for definiteness) . It can be shown along the 
lines of the proof of Lemma |3] that for all / G F, 

e(S, g (A»)J,e-^ 8 )<3e(S,A,f). 

Note that the cost is 0(u 2 ) (which is usually a logarithmic term) times the 
cost of computing the norm ||<?j — gj\\ (which depends on the structure and 
dimension of G, and on the - possible - sparsity of the tft, see also the 
discussion in Section 3 of jllj). 

Corollary 2. Let D,K, F C F(D,K), G, k G N and S,Si : F -> G 
(I = 0, . . . , k) be as in Lemma\^ Assume . . . , G N satisfy 



- 4 

1=0 

Let ng, . . . ,i%k G Nq and put n = Yli=o u i n i- Then 

k 

eUS,F)<2J2$ l (S l ,F). 

1=0 

If G = loo (T) , then the relation holds with constant 1 . 

This is an obvious consequence of Lemma|2]and Corollary^ In the sequel 
we need the following mappings. Let m* G N and define (3 : R — > Z[0, 2 m ) 
for z G R by 

[ if z < -2 m */2-i 

(3(z) = \ [2 m *l 2 {z + 2 m */2-l)j if _ 2 m*/2-l < z < 2 m */2-l ( 7 ) 

( 2 m * - 1 if z> 2 m *' 2 - 1 . 

Furthermore, let 7 : Z[0, 2 m *) — >■ R be defined for y G Z[0, 2 m *) as 

7 (y) = 2- m */2 2/ _ 2 ™*/2-l_ (g) 

It follows that for _2 m */2-i < z < 2^/2-^ 

7(/?W)<^<7(/3W) + 2- m * /2 . (9) 
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Proposition 1. Let D be a nonempty set and let / F C X C JP"(D,R) ; 
where X is a linear subspace equipped with a norm \\ \\x, such that 

(i) supj g p \f{t)\ < oo for each t £ D, and 

(ii) X separates the points of D in the following sense: Given to £ D and 
a finite subset D C D\{t }, there is an g £ X with g(to) / and g{t) = 
for all t G Dq. 

Let J : F — > X be the embedding map, let G be a normed space and 
S : X — > G a bounded linear operator. Then for all h, n G N, < 6±, 9 2 < 1, 

4+2n(SJ, F, 9 1 + 9 2 - 9 1 9 2 ) < e~(J, F, 6,) e«(S, B(X), 9 2 ). 

Proof. Let 5 > 0, let A be a quantum algorithm from F to X with q(A) < h 
and 

e(J, A, F,0!)< 63(7,^,0!) + $:= ffl . (10) 

Put 

a = (7i + (5. (11) 
Let A be a quantum algorithm from £>(X) to G with (7(A) < n and 

e(5,A,B(X),e 2 )<eS(5,B(X),e 2 ) + <J:=(72. (12) 

Let 

A = ((A l )*- \(b l )*- \lp), 

with 

4 = (Ql,(tf«)?=4>)> 

and 

Qi = (mi,fhi,m",Z h T h Pi). 
for Z = 0, . . . , k — 1. Furthermore, let 

^=((Mto 1 ,(&0to^), 

and for Z = 0, . . . , k — 1, 

A = (^,(^)Jio)> 

and 

Qi = {mi,m' l ,m",Z l ,Ti,(3 l ). 
We need some auxiliary functions and relations. Let 

D A = {n(i)\l = 0,...,k-l,i€Z l } 
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(the set of all points at which the quantum algorithm A queries the function) . 
By assumption (ii), for each t G Da there is a gt £l such that g(t) = 1 and 
g(s) = for all s G Dj±\{t}. Let Mi = ma,x t< z£) A \\gt\\x- By asumption (i), 

M 2 := max I f(t)\ < oo. 

f€F,t€D A 

Now choose the m* in the definition of the mappings f3, 7 in and ((HJ) in 
such a way that for a G R with |a| < M2, 

|a- 7 o/3(a)| < M^ 1 \D A \- 1 6. 
Define for / G F and x = (x , x~ k _ x ) G nfco Z [°> 2 ™0> 

= ( / - ^) + E (t /(*) - ft))* I 

\ t€D A J 

(recall that (f(x) G X). Then hf x G X, 

h Lx (s)=a-\ 1 o(3of( s )-{p(x)(s)) (s G Da), (13) 

and 

||/ - <p(x) - crh f>x \\x < Mil^lMf 1 !^!-^ = 5. (14) 

Moreover, 

||fy>IU = v^Mf ~ &( x )) ~ (f ~ &( x ) ~ ° h f,x)\\x 

< a-WKf-^Wx + d). (15) 

We build an algorithm as follows. It has k + k cycles. The first k cycles are 
exactly those of A. After the k measurements we have the result, say 

X (^0) x li • • • ) ^fc— 1)) 

(from which tp(xo, • • • > x t-r) would be computed - but we don't do that yet). 
Next the k cycles of A follow, with certain modifications. In each cycle we 
add rh = Y2i=o ™l <l UD its which are initialized in the state 

\x) = \x ) \xi) . . . |a?g_ 1 ) 

and remain there all the way. We add m* further auxiliary qubits, initially 
set to zero (and being zero again at the end of each cycle). We also want 
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to modify the queries Qi of A. For < / < k introduce the following new 
query: 

Qi = (mi + m + rn*,m'i,m*,Zi,ri,l3), 
where (3 was defined in ((7J). Define a unitary operator V\ on 

Hm't ® H m n <g) i/ m; _ m J_ m » (8) <g> il m * 

by setting for 

\i) \z) | it) |ar) |v) G C m j (g> C m » <g> C m ,_ m /_ m » (8) C„ (g> C m » 

Vi \i) \z) \u) \x) \v) = \i) \z j 9,((7- 1 (7(t;) - fa (*))))) |u) k) |i>> 
if i G Z[, and 

F { |i)k>l«>k>b) = |i)k)|«)|x)|«) 

otherwise (recall that (p(x) £ X C F(D, R) and r^(i) G D, so the respective 
expression above is well-defined). We also need 

Wi \%) \z) \u) \x) \v) = \i) \z) \u) \x) |e«) 

with = (2 m — v) mod 2 m . Now we consider the following composition 

PiQijPiWMPiQuPu (16) 

where P; exchanges the \z) with the \v) component. Let us look how the 
combination (|16|) acts, when the last m* qubits are in the state |0). Assume 
i G Zi. Then 

\i) \z) \u) \x) |0> 

is mapped by PiQijPi to 

|i)|*>|«>|*)|/?0/07Ki)>. 

Using (fT3*|) . we see that V\ produces 

\%)\z@ 0i (tr- 1 ( 7 o /3 o / o 7](i) - £(z)fa(i)))) ) |«> \x) \0 o f o tj(»)> 

= |i) |z e A ° \u) \x) |/3 o / o 7i(i)) . 

Finally, the application of PiQijPiWi leads to 

\i)\z®(3 l oh f]X o r l {{))\u)\x)\Q) = (Q l>hf Ji) \z) \u)) \x) |0) . (17) 
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It can be checked analogously that if % Z\, we also end in the state given by 
the right-hand side of (|17|). That means, the combination (|16j) acts as if we 
apply the original query Qi, but with / replaced by hf )X . Now we replace 
each occurrence of Qij by this string (fTo]) . while the Uij are replaced by 
Uij, which are the Uij, extended to H m (g) ® il m * by tensoring with the 
identity on H,^ <g> H m * . The respective 6/ are defined in such a way that 

bi(x , ■ ■ -,x- k _ v (y o ,x,0),.. . , (yi-x,x,0)) = (k(y , . . .,yi-x),x,0) 

for all x ,.. -,x^_ v y , . . . ,yi-i, l = 0,...,k-l. 

After the completion of the k + k cycles, let the measurement results be 

xq, ■ ■ (Vo,x,0), (y fc _i,x,0), 

where, as before, x = (xo, . . . , xt^)- Then we apply the mapping (p defined 
by 

tp(x, (y ,x,0), . . . , (y k -x,x,0)) := Stp(x) + aip(yo, . . . , Vk-x)- 
Denote the resulting quantum algorithm from F to G by B. Clearly, 

q(B) = n + 2n. (18) 

By (|17|). the modified A-part, applied to / G F, acts like algorithm A, 
applied to More precisely, in algorithm B, applied to /, given x as 

the outcome of the measurements of the first part of B, the probability of 
measuring 

(2/o,a5,0),... ,(y fc _i,x,0) 
in the second part of B is the same as that of measuring 

2/0, • • • ,2/fc-i 

in algorithm A, applied to For a fixed f £ F, we have, by (|10|) . with 
probability at least 1 — 0x, 

\\Jf ~ fi(%)\\x < o-l, 
thus, by (|T5|) (recalling also Jf = f), 

\\hf, x \\x < ^{oi + 5) = 1. 
Thus, for fixed / G F, with probability at least 1 — Ox, 

h f , x G fi(X). (19) 
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But for each x satisfying (fTU|) . we have by (|12|). with probability at least 

1-02, 

\\Sh fjX - ip(y ,... ,y k -i)\\ G < cr 2 . (20) 

Summarizing, we see that (jl9|) and (|2Uj) together hold with probability at 
least (1 - 6>i)(l - 6 2 ). We have 



< 
< 



\\SJf - <p(x, (y ,x,0), ... , (y k -i,x,0))\\ G 
\\SJf-S<p(x) - <T<p(y ,...,y k -i)\\G 
\\SJf - Scp(x) - crSh f)X + aSh f)X - a<p(yo, . . . ,y k -i)\\G 
\\S{f - <p(x) ~ oh fjX )\\ G + <j\\Sh ftX - tp(y , . . . ,yk-i)\\G 
\\S\\8 + aa 2 , 



by (dU) and (gDJ), with probability at least (l-0i)(l-0 2 ). Thus, using ifTTji. 

e(SJ,B,F,e 1 + e 2 -9 1 9 2 ) < \\S\\5+(a 1 + S)a 2 , 
hence, by (fT5|). lfTU|l. and JUJ), 



'n+2n 



< ||5||«5 + (4(J, F, + 2<y)(<£ (5, 2 ) + J). 

Since 5 > was arbitrary, the result follows. 
Corollary 3. Let u\,u 2 € N imtt 

e -i/i/8 _|_ e -v 2 /8 _ e ~(u 1 +u 2 )/8 < 

T/ien under the same assumptions as in Proposition^ 



□ 



(5J,F)<4 e ^(J, J P)e^(S,0(X)). 



Proof. By Proposition ^ and Corollary^ 



'i/in+2^2f 



n (5J,F) < e ^-(J,F,e-^/ 8 )e^ n (5,^(X),e^/ 8 ) 



< 4e3(J,F)eS(5,B(X)). 



□ 
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4 Approximation of Finite Dimensional Embed- 
dings 

For N G N and 1 < p < oo, let denote the space of all functions 
/ : Z[0, N) — ► R, equipped with the norm 



r JV 



1 

AT £ 



i=0 



if p < oo and 



rjv = max |/(i)|. 
0<i<N-i 1 v " 



Define Jp q : — ► L q to be the identity operator J^ q f = f (/ G ). 

~iN,M . r JV . r TV 



Furthermore, for a real M > define the operator C P q M : — > for 



/ = (/(O)S 1 as 



l^pg o otherwise. 



Lemma 5. Lei 1 < p, q < oo. There is a constant c > sucht that for all 
n, N G N, onrf M £ R imi/j M > 0, 

ea(qS^,B(^)) = 

lo/ieneuer 

M > c(iV/n) 2 / p max(log(n/\/iV),l) 2 / p . 

Remark. Throughout the paper log means log 2 . Furthermore, we often use 
the same symbol c, c\ , . . . for possibly different positive constants (also when 
they appear in a sequence of relations). These constants are either absolute 
or may depend only on p and q - in all statements of lemmas, propositions, 
etc. this is precisely described anyway by the order of the quantifiers. 

Proof. This is an immediate consequence of the proof of Proposition 1 and 
Corollary 3 of Namely, it contains an algorithm with n queries that 
produces, with probability > 3/4, all indices i with \f(i)\ > M and for each 
such i an (arbitrarily precise) approximation yi to fit), where M is any 
number satisfying 

M > c{N/nf/' p max(log(n/v / iV), l) 2/p . 

□ 
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Lemma 6. Let 1 < p < q < oo. For all N G N, and M G R toii/i M > 0, 



sup lu ^ M 



ll/-<' M /ll^<M 1 -^. 

Proo/. We have for / G B(Lp) 
1 ^ 



iV 

1 w 

< 



1 ^ 

- max i/(i) - (c»> M f)<j)r p E 1/(0 - (qS ,Af /)(oi p 



iV j 

M<? _p AT-1 
i=0 

□ 

Next we give an upper bound. 

Proposition 2. Lei 1 < p, q < oo. In the case p < q there is a constant 
c > such that for all n, N G N 

eS( J£, B(i# )) < c min I (- log (n/ViV + 2) J , iV 1 ^ 1 /* j . 

In £/ie case p > q, we have 

el{J^B{L»))<l. 

Proof. For p < q the estimate involving the first term of the minimum 
follows from the previous two lemmas, since by Lemma 6 (i) of [H], 

e«(J»,B(L»))<e«(C% M ,B(LZ))+ sup \\f - C% M f\\ L „. 

The estimate involving the second term is a trivial consequence of \\Jp q \\ = 
AT 1 /?- 1 /?. The case p > q follows from || = 1. □ 

Before we derive lower bounds we recall some tools from [B]. Let D and 
K be nonempty sets, let L G N, and let to each u = {uq, . . . , liL-i) £ {0, 1} L 
an / u G ^(-D, K) be assigned such that the following is satisfied: 
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Condition (I): For each t G D there is an I, < / < L — 1, such that 
f u (t) depends only on ui, in other words, for u,u' G {0, 1} L , U[ = u\ implies 

/«(*) = /«'(*)■ 

For -u G {0, 1} L let \u\ denote the number of l's in u. Define the function 
g(L, I, V) for L G N, < I + V < L by 



The following was proved in jHj, using the polynomial method [3] and based 
on a result from |16j : 

Lemma 7. There is a constant cq > such that the following holds: Let 
D, K be nonempty sets, let F C T[D, K) be a set of functions, G a normed 
space, S : F — > G a mapping, and L G N. Suppose (fu) u e{o,i} L — F(D,K) 
is a system of functions satisfying condition (I). Let finally < I ^ I' < L 
and assume that 

fu&F whenever |u| G {Z, I'}. (22) 

Then 

el(S, F) > \ min { \\S(f u ) - S(f u >)\\ \ \u\ = I, \u'\ = I'} (23) 
for all n with 

n< c q(L,1,1'). (24) 

For the case q = oo we give the following lower bound. 

Proposition 3. Let 1 < p < oo. There are constants c\,C2 > such that 
for all n,N G N with n < c\N 

e n( J p,ooiB{ L p )) > C2 min ^y P ,N^ . 
Proof. First we assume 

n < c \Hv, (25) 
where cq is the constant from Lemma Put 

L = N, 1 = 0, I' = 1. 

Then _ 

n<coVL = cog(L,l,l'). (26) 
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N l/p if i = j 

otherwise. 



Define ipj (j = 0, . . . , L — 1) by 

if>j(i) = 

Note that tf)j G B{L%) and 

II^oo^-IIl- =H^llL-=iV 1/p . (27) 
For each u = (uq, . . . , Ul-x) G {0, 1} L define 



L-l 



f u = J2^r (28) 
i=o 

Since the functions V>j have disjoint supports, the system (/u)ue{o,i} i satis- 
fies condition (I) . Lemma [7| and relations and (|2*T)) give 

eliJ^B^)) > ^mm{\\J^f u -J^f ul \\ L N\\u\ = 0,\u'\ = l} 

= -N l lv. 
2 

This proves the statement in the first case. Let 

Cl = c /VV2. (29) 

Now we assume 

c Vn <n< ciN. (30) 

We set 

L = N, 1= \2cq 2 n 2 N- 1 ] , l' = 1 + 1. (31) 
It follows from j3UJ) that I > 2. Moreover, from IpTTJl. 



n < co^lN/2 (32) 
and, taking into account that I > 2, 

J/2 < I- 1< 2 C Q 2 n 2 iV" 1 , 

hence, by (JUJ) and pTjl. 

Z + 1 < 3Z/2 < 6cq ViV^ 1 < 6c 2 c 2 iV = JV/2. (33) 
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We have, bv l|33)> . (j33| ) and l|3T|l . 



n < c 0V W/2 < c min ^ j{N - j) < c g(L,l,l'). (34) 
j=l,l+l 



Now we define tpj G (j = 0, . . . , L — 1 



as 



Then 



(I + l)- l /PN x /v if i = j 
otherwise. 



= = {l + ^r 1/p N llp . (35) 

Defining the system (f u ) u e{o,i} L as above, it satisfies condition (I), and 
f u G B{Lp) whenever |u| =1,1 + 1. Lemma0 relations (jBU), (j33|) . and the 
left and middle part of (|3*3)) give 



1 

2 

i(/ + IJ-Vpjv-Vp > I(6cq 2 n 2 iV- 1 )- 1 / p iV 1 / p 

2/p 

^0 n -2/ PAr 2/ P _ 

2 • 6Vp 



□ 



The previous results together with Corollary El give lower bounds also 
for arbitrary q. 

Proposition 4. Let 1 < p,q < oo. There are constants cq, c\ > such that 
for all n, N G N, n < coiV i/ie following hold: If ' p < q, then 



e^,£(^))> Cl min( (3) 



N -,2/p-2/q 

n J 



log (n/ViV + 2)) 2/ ",iV 1 /p- 1 /^ ; 



and if p > q, then 



el(J"B{L"))> Cl {\ogN + ir^. 



'pqi 

Proof. Fix any 1/1 , 1/2 G N with 

e -i/i/8 _|_ e -^2/8 _ e -(i>i+i/ 2 )/8 < -w^ 

By Corollary El 

e (, 1+2 , 2 )n(^oo^(^)) < 4 e ^(J^,^))e^(J 9 Ar 00 ,e(Lf )), 
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therefore, 

e*(J»,B(L»)) > A-'el^JJ^BiL^eUJ^B^ ))~\ (36) 
By Proposition |2 

eUJ^oo, B(I%)) < c min ^ ^ log (n/ViV + 2) ) ^ , iV 1 /^ , (37) 

while by Proposition for n such that (y\ + 2z^)n < c\N {c\ the constant 
from Proposition |3J) 

^w-s..^)) * ^ min (((^k^) 2/P ^ 1/p ) 

> cmin ((—) /P ,NM . (38) 



We first consider the case n < min y/N,c\N/(v\ + 2^2) j- Then (|36j) . 1)37(1 . 
and (|3*%)) give 

e n(JpqiB(Lp )) > cN 1 /' p ~ 1 / q . (39) 
In the case yN < n < c\N/(yi + 2^2) we obtain similarly, 

e^,S(L^))> C (jIJ ^_iog(n/v^V + 2)J . (40) 

For p < q relations ()39() and (|4L)|) give the lower bound. In the case p > q 
we note that it suffices to prove the lower bound for n = \c\N/(vi + 2z^)J. 
But in this case (|40|) gives 

e^(J^,i3(^))>c(logiV + l)- 2 ^. 

□ 

To present the results in a form which emphasizes the main, polyno- 
mial parts of the estimates and suppresses the logarithmic factors, we in- 
troduce the following notation: For functions a, b : N 2 — > [0, 00) we write 
a(n, N) Xi og b(n, N) if there are constants Co, c±, c-i > 0, no S N, «i, 02 G R 
such that 

ci(log(iV + n)) ai b(n,N) < a(n, N) < c 2 (log(iV + n)) a2 b(n, N) 

for all n, N 6 N with uq < n < CoN. In this notation, we summarize the 
results of Propositions El and in 
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Theorem 1. Let 1 <p,q< oo. If p < q then 

( / M \2/p-2/q \ 

e*(J»,B{L»)) x log min f (-J , JV^M , 

and if p> q, then 

e n(.Jpqi B(L p )) ><log 1. 

Although the polynomial part of the order has been determined, there 
remains some room for improvements of the logarithmic factors. In the 
sequel we present some improvements of the lower bounds for particular 
situations. They involve a different way of applying Corollary |3J which is 
interesting in itself: we combine it with known results about summation. 
Furthermore, we will use these bounds for the Sobolev case We need 
to recall some previous results about summation. Define Sn : — > R by 

1 

i=0 

Let us recall what is known about the minimal query errors of the Sn- 

Proposition 5. Let 1 < p < oo. There are constants cq,ci,C2 > such 
that for all n, N G N with 2 < n < c N, 

cm' 1 <el(S N ,B(L*))<c 2 n- 1 if 2 < p < oo, (41) 
cm" 1 < e*{S N ,B{L%)) < c^rT 1 log 3 / 2 n log log n, (42) 

and 

cimin^-^-^n-^iV 2 /^ 1 ) < e%{S N , B{L*)) 

< c 2 min^- 2 ^- 1 ^, n- 2 / p iV 2 / ?, - 1 )(log(n/v / iV + 2)) 2 / p ~ 1 (43) 

ifl<p<2. 

The upper bound in the case p = oo is contained in 0, [2], the lower 
bound for p = oo is from |16j . while (|41|) for 2 < p < oo and 1)42(1 were 
obtained in [Hj, and (|4.3|) is from Jl]. For the case p = 2 we need another 
upper estimate, which is more precise than (|42j) for n close to N. It can be 
found in |TQ|, Lemma 6. 
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Lemma 8. There is a constant c > such that for all n, N £ N with n < N , 
e*(S N ,B(L%)) < cn-H(n,N) 3 / 2 logl(n,N), 

where 

l(n, N) = \og(N/n) + log log(n + 1) + 2. 

Proposition 6. Let 1 < p, q < oo. There are constants cq, c\ > such that 
for all n, N G N with N > 4 and n < c iV 

el{J» v B{L%)) > Cl */ 2<g<oo 
el{J^B(L%)) > dOoglog^-^aogloglogiV)- 1 
e%(J» q ,B(L»)) > cxOogiV)- 2 /^ 1 if \<q<2. 

Proof. The proof is similar to that of Proposition Fix v\ , V2 £ N with 



e 



-1/1/8 



_|_ e -^2/8 _ e -(n+i^)/8 < 



Since is the identity, we have SV = SpjJpq, and we deduce from Corollary 

El 

e (, 1+2 , 2 )n(^(^)) < *<%(j£,B(L?))t*(S N ,B{I<?)) 
which gives 

($(J* > 4~ 1 e^ i+2iy2)n (^^(L^))e^(^,e(Lf 

It remains to apply PropositionEJand Lemma[E]with n = [co(^i + 2i>2)~ 1 N\ . 

□ 

5 Comments 

First we discuss the cost of the presented (optimal with respect to the num- 
ber of queries) algorithms in the bit model of computation. For this purpose 
we assume that n and N are powers of 2. The algorithm behind Proposition 

n 



121 for approximating Jil is nontrivial only for n > yN . In this case it uses 



pi 

quantum queries to finds all coordinates of / 6 B{Lzf) with absolute value 
not smaller than 

M = c(iV/n) 2/p max(log(n/v / iV), l) 2/p 

(with c a concrete constant, see |14j). It follows from H/H^jv < 1 that there 
are at most 

NM~ P = 0{n 2 N~ l maxOog^/v 7 ^), 1)~ 2 ) 
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of them. The algorithm of finding them, which is described and analyzed in 
[Tl| . needs O(nlogiV) quantum gates, O(logiV) qubits, 

0{n 2 N~ l max(log (n/ViV), 

measurements, and 0(n 2 N~ 1 log N) classical bit operations. So the total bit 
cost is O(nlogiV) = O(nlogn), hence, up to a logarithm, the same amount 
as the number of queries. 

In the following table we summarize the results of this paper on the 
quantum approximation of the identical embeddings J^ q : B(Lp) — > 
and compare them with the respective known quantities in the classical 
deterministic and randomized settings. We refer to |7j and the bibliography 
therein for more information on the classical settings. The respective entries 
of the table give the minimal errors, constants and logarithmic factors are 
suppressed. We always assume n < cN. 











deterministic 


random 


quantum 


l <p < q< 


oo, 


n < a 


/N 


jyl/p-l/q 


jyl/p-l/g 


jyl/p-l/q 


l <p < q< 


oo. 


n > 


v/iV 


jyl/p-l/g 


jyl/p-l/g 




1 < q < p < 


oo 






1 


1 


1 



We see that the quantum rate can improve the classical deterministic 
and randomized rates by a factor of order iV -1 (for p = 1, q = oo, and n 
of the order of N). It is this case which will lead to a speedup for Sobolev 
embeddings by an exponent 1, see We observe that there are also 

regions where the speedup is smaller or there is no speedup at all. 
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